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ABSTRACT 

Aims. To demonstrate the time-approach to equilibrium of Ha-formation and protonation in models of diffuse or H I interstellar gas 
clouds previously published by the author. 

Methods. The microscopic equations of H2 -formation and protonation are integrated numerically over time in such a manner that the 
overall structures evolve self-consistently under benign conditions. 

Results. The equilibrium H2 formation timescale in an H I cloud with N(H) « 4 x 10^" cm"- is 1 - 3 x 10' yr, nearly independent of 
the assumed density or H2 formation rate on grains, etc. Attempts to speed up the evolution of the H2-fraction would require densities 
well beyond the range usually considered typical of diffuse gas. The calculations suggest that, under benign, quiescent conditions, H2 
in the diffuse ISM formation of H2 is favored in larger regions having moderate density, consistent with the rather high mean kinetic 
temperatures measured in H2, 70-80 K. 

Formation of is essentially complete when H2-formation equilibrates but the final abundance of H^* appears more nearly at the 
very last instant. Chemistry in a weakly-molecular gas has particular properties so that the abundance patterns change appreciably as 
gas becomes more fully molecular, either in model sequences or with time in a single model. One manifestation of this is that the 
predicted abundance of is much more weakly dependent on the cosmic-ray ionization rate when n(H2)/n(H) S 0.05. In general, 
high abundances of do not enhance the abundances of other species {e.g. HCO^) but late-time OH formation proceeds most 
vigourously in more diffuse regions having modest density, extinction and H2 fraction and somewhat higher fractional ionization, 
suggesting that atypically high OH/H2 abundance ratios might be found optically in diffuse clouds having modest extinction. 

Key words, interstellar medium - molecules 



1. Introduction. 

The fact that individual H I/diffuse clouds have a substantial 
component of molecular hydrogen has been recogniz ed observa- 
tionally since the original Copernicus observations ( Jurai 119741: 
ISpitzeit Il978l) and a quite high fraction of hydrogen in the 
nearby di ffuse interstellar medium (ISM) overall is in molecu- 
lar form jSavage et all Il977t iLiszt & Lucasl |2002|) . More sur- 
prising perhaps is the recent discovery that this molecular 
component hosts a ric h polyatomic chem istry with relatively 
large amounts of jMcCall et al.L 200Z) and a dozen o ther 
species including such ioii s as HCO+ (lLucas&Lisztl[T996h and 
HOC^ d Liszt et al.', '2004) and molecules as com plex as 3h2 
jCox et a l., 1988; Lucas & Liszt, 2000), H2CO (L iszt & Lucasl 
Il995t iMoore & Marscheit 119951: ILiszt et all {20^ and NH3 



(lNashl[T990; Liszt et al.Ll2002h . 

The presence of copious H2 and relatively large amounts 
of in H I clouds is supported theoretically when self- 
shielding (iLee et al.L 19961) is inc l uded i n equi l ibrium mod- 
els of diffuse gas dLiszt & Lucasl l2000l [Lisztl l2003l) using 
empirically-determined H2 formation rates (lJuraL ll974t l ISpitzeii 
[1978; Gryetal., 2002); observed HD/H2 ratios are also ex- 
plained as long as the cosmic-ray ionization rate is taken large 
enough to support the inferred proton density in the presence of 
atomic-ion neutralization on small grains. However, formation 
of H2 is a notoriously slow and supposedly rather fragile pro- 
cess. Moreover, it remains to be determined whether any such 
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equilibrium is actually attained in the interstellar medium, which 
is the subject of this work. 

Here we show the approach to equilibrium in the models dis- 
cussed by Liszt & Lucas (2000), Liszt (20dl) and|Liszt (2003*): 
we present calculations of time-dependent H2 and Hs^ forma- 
tion in small, mostly-static gas parcels meant to represent typical 
diffuse "cloud s", like, for instance, the "standard" H I cloud of 
ISpitzeii ([T978h . As we discuss, the distinguishing characteristic 
of the chemistry in such clouds is that the H2 formation process 
is unsaturated (n(H2) < n(H)/2). For this reason the equilibrium 
timescales for H2 formation cannot be calculated from simple 
scaling arguments based on the local microscopic H2 formation 
rate, nor are they shortened by self-shielding, or in many cases 
by assuming higher density or even a higher rate constant for 
H2 formation on grains, etc. Furthermore the abundance of H3^ 
does not always vary as might be expected from consideration 
of its equilbrium chemistry in more fully molecular gas. 

Section 2 gives details of the methods used in the calcula- 
tions, results of which are presented and discussed in Sect. 3 and 
4. 

2. Details of the calculations and microphysics 

2.1. Framework and initial conditions 

The calculations presented h ere build on a frarn ework which 
was established in the work o f lWolfireet an(ll995l) . Specifically, 
the two-phase model of heating and cooling developed by 
[Wolfireetal. (1995) (see also Wolfire et al. (2003)) was ap- 
plied to small, uniformly illuminated spherical gas clots of con- 
stant density, intended to represent diffuse 'clouds'; our model. 
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howev er, uses the gas-phase abundances of ISavage & Sembaciil 
( 1199 6') which makes them a bit hotter. Assumption of a geometry 
is required because the accumulation of H2 and CO depends on 
shielding by dust and other molecules and is, therefor e, strongly 
non-li near. We used the shielding factors calculated bv lLee et al.l 
(Il996l) . Typical densities and column densities considered here, 
n(H) - 32 cm"^, N(H) = 4 x 10 ^° cm"^ correspond to a Spitzer 
"standard cloud" ( SpitzeE [l928l) having a dimension of a few pc 
and median reddening 0.05 mag over the face of the cloud. The 
central column density in such an H 1 cloud corresponds to the 
regime where the appearance of high Ha-fractions is first noted 
in the ISM. 

In general, the hydrogen densities in our models (number or 
column) are assumed values. The temperature and thermal pres- 
sure follow from the details of the calculation, chiefly through 
the assumed gas-phase abundances of the dominant coolants 
carbon and oxygen, and the influence of the radiation field in 
heating the gas via the photoelectric effect on the same small 
grains which are so largely res ponsible for th e overa ll charge 
balance. As shown in Fig. 9 of iLiszt & LucasI (l200 0b, a given 
value of the scaling factor (Go) for the overall radiation field 
produces a roughly constant thermal pressure in the outer re- 
gions of our models, increasing somewhat with assumed density 
and decreasing slightly going inward from the cloud edge. The 
generally-accepted value for the thermal pressure of the ISM, p/k 
2 - 3 X 10^ cm-^ K (Jenkins, 2002; Jenkins & Tripp, 2001) in- 
deed applies to our typical model of the Spitzer "standard cloud" 
and the very densest models considered here have pressures a 
factor 2-3 times higher, as seems observationally to be the case 
as well. 

The initial condition here is simply the equilibrium diffuse 
cloud model which would be calculated in the absence of any 
molecular chemistry, as is more ty pical in discussions of the 
overall structure of the diff'use ISM ( Wolfire et all l2003h . That 
is, the initial model is equilibrated and self-consistent in terms 
of thermal and ionization equilibrium, but no molecule forma- 
tion has occurred. Our models include the rudiments of molec- 
ular chemistry as discussed in several recent papers concerning 
question s of equilibrium CO, H2, HP and Hj^ f ormation in dif- 
fuse gas dLiszt & Lucasll2000l:IUsz3.l2002ll2003 ). Incorporation 
of the chemistry should not normally upset the basic microphys- 
ical balance (ionization, thermodynamic, etc.) in diffuse gas be- 
cause such small fractions of oxygen and carbon are sequestered, 
although inclusion of OH and H2 formation can affect some as- 
pects of the distribution of charge among species. 

Another effect mitigating the influence of H2 formation in 
diffuse gas is the dominance of atomic-ion neutralization by 
small grains, substantially lowering the and electron den- 
sity which would otherw ise obtain in the presence of cosrn ic-ray 
ionization of hydrogen dWolfire et all 119951; iLiszj |2003|) . This 
effect is of the utmost importance in understanding the inter- 
play between the cosmic ray ionization rate (h ere taken as 
the rate per free H-nucleus) and abundance of Hj^ ( Liszt, 2003^) 
. Detailed calculation including grai n charging (Draine & Sutin, 
119871; IWeingartner & Drainel 1200 lb is required if the heating 
ra te is to be properly cal culated according to the prescription 
of lBakes & Tielensid 19941) . 

2.2. Formation of H2 

In this work we distinguish between the fonnation of H2, which 
we consider to be a microscopic physical effect occuring on indi- 
vidual large grains and involving individual or paired H-atoms, 
and the accumulation of H2 in the ISM, which is a macroscopic 



effect involving assumptions about cloud geometry, ISM topol- 
ogy, and so forth. The purpose of this work is to understand the 
macroscopic aspects of the accumulation of H2 in the ISM and 
we employ only the most standard, empirical, microscopic de- 
scription of its formation. 

That is, the volume formation rate of H2 is expressed as 
dn(H2)/ dt ^ n(H ) n(H I ) Rq <V7->= 3 x lO '^cm^ s ' n(H) n(H 
I) ^/T^ dSpitzeii Il978l) where flie factor n(H) (the total density 
of H-nuclei) represents a constant number of large H2-forming 
grains per H-nucleus in the gas; n(H I) is the density of atomic 
hydrogen (after protons, H" ions, and H2, etc. are reckoned) and 
the thermal velocity and/or kinetic temperature terms correspond 
to the speed at which H-atoms move through the gas, possibly 
encountering large grains. The rate constant Rg nominally in- 
cludes a sticking coefficient. 

The inferred rate in the diffu se ISM has been r e mark- 
ably stead y over the last 25 years dJural 119741; ISpitzeii 119781; 
iGrv et al.L l2002h and use of this simple formulation, com- 
bined with the modern shielding coefficients, accounts ex- 
tremely well fo r the observations of H2 in diffuse lines of sight 
dLiszt & Lucasll20b 0). As discussed in Sect. 3 (see Fig. 4) vary- 
ing the assumed value of Rq by factors of a few has appreciable 
effect upon the amount of H2 which is produced but the equilib- 
rium time scale changes little. As an aside we note that the gas- 
phase formation of H2 via H , which could account for some of 
the very low molecular fractions seen along lines of sight having 
low extinction (Liszt, 2002), is included in the present calcula- 
tions as well but without much effect. 



2.3. Solving the H2-accumulation problem 

The equilibrium solutions for the radial distributions of CO, H2 
and H3^ described in our earlier papers were achieved through 
what is essentially a relaxation method. An initially atomic 
sphere is divided into a substantial number (typically 100) of 
thin radial shells and the molecular densities are calculated 
for each, working inward, self-consistently accounting for the 
shielding which accrues from each configuration. Each time a 
quasi-equilibrium state is reached for the entire structure (that is, 
when the calculation has successively converged in each shell), 
the calculation starts anew at the outermost shell (whose H2- 
abundance was, after all, calculated in at least partial ignorance 
of the conditions interior to it) and the process is repeated until 
the whole structure varies sufficiently little between iterations. 
The process is not particularly difficult or unstable and conver- 
gence is rapid, typically requiring only a few seconds on a recent 
laptop computer. Note that the thermal and charge balance, etc. 
are continually updated during the calculation. 

The time-dependent calculations basically just wrap the 
equilibrium solver in an outer loop which allows the entire struc- 
ture to evolve self-consistently over short times. The time step, 
which initially cannot be as large as even 1000 years in most 
cases, is examined and lengthened by a factor of V2 whenever 
the structure changes sufficiently little between time steps. The 
full time-dependent calculation is somewhat tedious, typically 
requiring several hours per model and much exercising of the 
tiny (but assertive) fan on the author's laptop computer in cases 
of high H2-fraction. 
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2.4. Elementary analytic considerations of the 
time-dependent solutions for H2 

Although the real solution to the time-dependent accumulation 
problem is strongly non-linear and can only be achieved globally 
over an entire geometrical construct, several important aspects 
of the problem may be understood with reference to a simple 
parametric treatment of purely local behaviour 

Accordingly, the local rate of accumulation of molecular hy- 
drogen may be written 



dn(U2)/dt = -n(H2)rH2 + ^Gn{H)n{HI) ^lT^ (1) 

where Fjj^ is the rate of H2 destruction (mostly by cosmic 
rays and ultra-violet photons, although interaction of He^ and 
H2 is included in the numerical calculation), Tk is the kinetic 
temperature (calculated in the models) and Rg = 3 x 10"'^ cm"* 
s^' (see Sect. 2.2). After rewriting n(H I) = n(H) - 2 n(H2) (for 
the present purposes; all forms of hydrogen are considered in the 
numerical models) the solution is 

«(H2)|, = bla + («(H2)|,=o - ^/fl)e-"' (2) 

where a = T^^ + ^^Gn(H) ^|T2, b = Rcn(H)^ a/Tk". 

The time constant in the problem is 1/a, where a is a lin- 
ear combination of the per-particle formation and destruction 
rates, observing the conservation of H-nuclei, and is therefore 
relatively short in weakly-shielded gas. The rate of accumula- 
tion of H2 takes on its highest value dn(H2)/dt = b at time and 
the equilibrium solution n(H2)|,^oo - b/a may be thought of as 
the result of accumulation at the maximum rate over a period 
of time 1/a. Thus, self-shielding (lowering Tj^^) cannot increase 
the rate at which H2 forms or speed up the accumulation of H2, 
it only increases the final amount of H2 by allowing the accu- 
mulation process to work longer. By lengthening the time-scale, 
self-shielding may also provide a buffer against sudden change 
once equilibrium is attained. 

Increasing the microscopic formation rate constant or den- 
sity will hasten the pace at which H2 accumulates but even so, 
global equilibria are not generally achieved more rapidly in cases 
where the H2 fraction would otherwise be small, as discussed in 
Sect. 3 (see Fig. 4); it only happens that more H2 is made (which 
requires somewhat more time). Maintenance of high equilibrium 
H2-fractions requires Fj^^ ^ f^Gn{H) VIk, which is roughly 
equivalent to reducing the photodestruction rate (of order 10 
s"' in free space) to insignificance i.e. Fjj^ ^ ^jj. 

2.5. Weak vs. slow processes 

Indeed, the situations which equilibrate most rapidly are those 
for which the H2 abundance is smallest, which exemplifies the 
fact that even processes like H2 formation or cosmic-ray ion- 
ization - which seem weak because they have low rates - may 
equilibrate quickly when not called upon to do very much; the 
time constant is then determined by the overall destruction rate, 
which may be high. It is for this reason that formation of Hs^ - a 
much higher-order process which after all requires formation of 
H2, cosmic ray ionization of H2 and further interaction of H2^ 
and H2 (see Sect. 4.2) - is essentially complete as soon as H2 
accumulation equilibrates: one need not wait for a further period 



2.6. The cosmic-ray ionization rate 

The evolution of the H2 fraction is largely independent of the 
assumed value of in our models but the standard cosmic-ray 
ionization rate assumed in the models in this work is - 10"'^ 
s"'(the cosmic-ray destruction rate per H2 is slightly more than 
2^//). This is approximately one order of magnitude larger than 
what is often used in calcu lations of cosmic-ray ionization of the 
ISM jWolfire et all 1 19951) especially for dense dark cloud chem- 
istry and is about that which is needed to explain the overall 
abundance of Hs^ (see Sect. 4). The equilibrium concen- 
tration is nearly proportional to in our models having larger 
H2 fractions (see Fig. 5) owing to the intervention of small-grain 
neutralization of protons: the formation can be pushed harder 
without unduly increasing the electron density and recombina- 
tion rate for H3^. However, for reasons discussed in Sect. 4, the 
concentration of is more weakly dependent on in regions 
of small H2-fraction or when is very large. 



2. 7. Heating by cosmic rays and H2 -formation 

The current assumption is that the diffuse ISM is heated by 
the photoelectric effect on small grains in both the warm phase 
(where x-rays make a larger contribution) and cool phase. The 
heating rate as a func tion of density in partly-shielded gas is 
sho wn in Fig. 3 of Wolfire et al.l (1 1995b based on the work of 
Bak es &Tielensl(ll994 . As may be seen from that Figure, the 
cosmic ray heating rate is roughly 1 .2 dex below that of the pho- 
toelectric effect in diffuse gas and 2-2.5 dex lower in cool gas, 
for ^//= 10^'^s"'. Thus the higher cosmic-ray ionization rates 
considered here may imply a substantial contribution to the heat- 
ing of the warm gas by cosmic rays, barkening back to the very 
earliest discussions of two-phase equili brium, before the role o f 
diffuse X-ray heating was incorporated dGoldsmith et aUri969l) . 
The highest cosmic -ray rates considered here still do not pro- 
duce dominant heating effects in the cool neutral gas, but much 
hi gher ionization rates 4"// ^ 10"'^s"' were recently suggested 
by iLe Petit et"an(l2004 " for the line of sight to ^ Persei. 

In our models there is a modest contribution to the heating 
due to H2 formation, comparable to that arising from heating by 
cosmic-ray ionization of H I. The volume heating rate due to 
the latter is n(H I) (Ecr) where (E„) ~ 8eV is the energy 
available for heating after the primary ionization. The equivalent 
rate for H2-formation is RGn(H)n(HI) ^/TK{Ell^) where (Ejj^) 
represents that fraction of the heat of formation of H2 which is 
available to heat the gas. The heating rates for the two processes 
are very nearly equal under typical conditions, = 10"'^ s"', 
Tk = 80 K, n(H) = 32 cm-^ and (E^^) = 1-5 eV. 



2.8. The recombination rate of H^"" and other chemical 
reactions 

The chemical rate constants used in our model s have been 
taken from the UMIST database dLeTeuff et all l2000b (see 
http://www.rate99.co.uk) including, for consistency with the re- 
sults of Liszt (2003), the rate constant for e -H recombination 
o(Tk) = 4 X I ^|T^ cm-' s"'. This rate has been measured 
m any times, inc l uding , most recently and probably accurately, 
bv lMcCalletaTI (120031) . Their measured values, at 23 K and 300 
K, are 70% of those used in our work. 
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Fig. 1. Density of molecular hydrogen n(H2) (left) and H^* (right) as functions of time within spherical, uniformly illuminated, 
initially atomic gas spheres of constant density of H-nuclei n(H) = 32 cm"-*, threaded by a cosmic ray flux yielding a primary 
ionization rate per H-atom = 10"'^ s"' . Three curves corresponding to shells at fractional radii 0, 1/3 and 2/3 from the center are 
shown for each of a series of models with differing central column density of hydrogen 0.5 x 10^" cm"^ < N(H) < 16 x 10^" cm"^. 
The physical size of the models in pc and their central (edge-to-edge) column density in units of 10^*^ cm"^ are equal and shown at 
the right of each curve. 



3. Time-dependent I-I2 accumulation in diffuse 
clouds 

3.1. Evolutionary effect of varying column density 

Figure 1 shows the evolution of the H2 and densities in a 
series of models, each having the same uniform total density of 
hydrogen n(H) - 32 cm"^ threaded by a cosmic ray flux yield- 
ing - 10 s '. The central column density of hydrogen 
across the clouds N(H) (from edge to edge) varies from 0.5 to 16 
xlO^'^ cm"^, as indicated at the right hand side of each curve. The 
physical size of the models (diameter) in pc is, coincidentally, 
equal to the column density as labelled in units of 10^" cm"^. 
Three curves are shown for each model, corresponding to the 
shells centered at fractional radii 0, 1/3 and 2/3 from the cen- 
ter. These curves separate at larger times when the outermost 
regions have reached equililbrium. For a sphere, the median line 
of sight occurs at a fractional impact parameter of 2/3 where 
the column density traversed is 3/4 that at the center. The me- 
dian line of sight through the model with central column density 
4 X 10^" cm ~^, 3 X 10^° cm" ^ corresponds exactly to the standard 
H I cloud of lSpitzed (Il978h . 

For the optically-thinnest model there is a large radial gra- 
dient in the H2 density and a wide disparity between the times 
required to reach equilibrium (10** - 10^ yr). The fact that such 



appreciable (easily detectable) amounts of H2 form in as little 
as 10'* yr demonstrates why it is so hard to find lines of sight 
having measurable reddening which are devoid of H2. For N(H) 
as large as 16 x 10^'^ cm"^ the H2 accumulation process is thor- 
oughly saturated (n(H2) n(H)/2) when equilibrium is reached 
at very large times approaching 10** yr. Models in which the 
equilibrium H2 -distribution is more uniform also evolve more 
uniformly, as illustrated in Fig. 2 which shows (at top) the time- 
dependent evolution of the radial variation of n(H2) in the model 
resembling a "standard" H I cloud. The H2 -density is remark- 
ably constant for fractional radii of 90% or more for most of the 
first 1 - 3 X lO*" yr, after which a more molecular core develops, 
as well as a very thin (actually, unresolved) transition layer into 
a more nearly atomic envelope. 

Clearly, self-shielding is important only at somewhat later 
times. In some models with moderate density, n(H2) behaves as 
a single power-law from center to edge, as illustrated in Fig. 2 
at bottom which shows the equilibrium solutions for n(H2) = 
32 cm"^ corresponding to the models used in Fig. 1. However 
this is not a general feature of the equilibrium solutions. Figure 
2 also illustrates the effect of finite resolution on the structure 
of the atomic envelope, which fractionally is substantially larger 
with higher resolution. The results in Fig. 1 are not affected by 
the illustrated doubling of the number of radial shells. 
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Fig. 2. Radial and secular variation of molecular hydrogen den- 
sity within models illustrated in Fig. 1. Top: Time evolution 
of n(H2) for n(H) = 32 cm-^ N(H) = 4 x lO^" cm^^. Labels 
at right give the elapsed time in units of 1000 yr Bottom: 
Radial variation in equilibrium for n(H) = 32 cm"^ and N(H) 
= 1,2,4, 8, 16 X 10^° cm~^, for calculations employing 96 and 
192 radial zones. 



As shown in Fig. 1, even a standard H I cloud model has 
a very substantial fraction of molecular hydrogen in its inte- 
rior. This of course is borne out observationally (lJurai 1 19741) as 
molecular fractions of 25% -45% are derived for even the dif- 
fuse ISM observed in Mv/optical absorptio n lines of H2 and CH 
dSavage et al.L 119771 : iLiszt & Lucasl |2002|) . Less obvious, how- 
ever, is that the conditions so favourable to H2 may actually work 
for as long as the calculations require to reach equilibrium in Fig. 
1, i.e. for times as long as 10^ yr. 



3.2. Evolutionary effect of varying number density 

It seems natural to inquire in what ways the evolution might be 
hastened. Shown in Fig. 3 are results for series of models of dif- 
ferent number density 8 cm""* < n(H) < 128 cm""* at column den- 
sities N(H) - 2 - 16 X 10^" cm"^. At least for moderate column 
densities, as shown at the left in Fig. 3, increasing the number 
density within the range expected for diffuse gas has surprisingly 
little effect on the time required for the models to reach equilib- 
rium: the increased rate of accumulation is compensated by the 
tendency to produce more H2, making the equilibrium timescale 
about constant over a wide range of density for the standard H I 
cloud column density. 

Only when very high H2 fractions are attained at a given den- 
sity does increasing the density shorten the equilibrium time (i.e. 
Fig. 3 at right, or for the two very highest densities at left). With 
reference to Fig. 3, note that, although the local accumulation 
rate initially varies as n(H)^ VTk^ n(H)^^^ in thermal pressure 
equilibrium, the total amount of H2 produced within a model at 
any time is much less dependent on density because the more 
tenuous models are larger, varying as l/n(H). Obviously, all 
models in which the H2 accumulation saturates at high molecu- 
lar fractions produce the same N(H2) if they have the same N(H), 
so less dense models having the same N(H2) then make more H2 
molecules overall because they are larger. 

3.3. Does size matter? 

Physically larger clouds tend to produce more H2 in both series 
of models just discussed. Those in Fig. 1 with larger column den- 
sity N(H) are physically larger and have higher internal molecu- 
lar number density, so, at fixed number density, physically larger 
regions of higher column density are responsible for the very 
great bulk of the accumulation of H2. In Fig. 3, molecular num- 
ber densities are higher in clouds of higher n(H), but usually not 
by such large amounts as to compensate for the very large dif- 
ferences in volume. That is, the model volumes vary by a factor 
(128/8)^ ~ 4000 whereas n(H2) varies by factors of 100-500 for 
the models with N(H) = 4-16x10^° cm"^. Over a substantial 
range in N(H), it is expected that the bulk of the molecules will 
form in (again) physically larger regions, but of lower density. 

So is the bulk of the interstellar H2 really formed in large ten- 
uous regions, or in more compact, denser ones? Observationally, 
the mean column density-weighted kinetic temperature inferred 
for H2-form ing diffuse gas is 70-80 K dSavage et al.l 119771 
IShull et aLl 12004.) . which, for typical interstellar thermal pres- 
sures, implies densities like those of the standard H I cloud Q. 
Attempts to form the observed H2 more rapidly in regions of 
high density must accou nt for these high kinetic temperatures; 
perhaps it is relevant that lJenkins & Trippl(l200ll) found a small 
but pervasive fraction of material at much higher than normal 
thermal pressure in typical diffuse clouds probed in the fine- 
structure excitation of C I. This high-pressure component is gen- 
erally ignored in calculation of the mean thermal pressure in the 
ISM, but would approximately double it. 



3.4. Changes in formation and dissociation rate constants 

The empirically-inferred microscopic formation rate of molec- 
ular hydrogen formation has chang ed r e markably little since 
H2 was first observed in the ISM dJural 119741: ISpitzeA 119781: 



' Note that the surveys from which these mean temperatures are de- 
rived have somewhat low sample mean H2-fractions, 0.17-0.18. 
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Fig. 3. Density of molecular hydrogen n(H2) as a function of time within spherical, uniformly illuminated, initially atomic gas 
spheres of differing density of H-nuclei n(H) = 8-128 cm"^ (as shown at the left of each family of curves), threaded by a cosmic 
ray flux yielding a primary ionization rate per H-atom (h - 10 s '. At right, all models have a central column density N(H) 
= 16 X 10^° cm"^ (from edge to edge); at left, results are interspersed for N(H) - 2 (dashed; red in appropriate media) and 4 
xlO^° cm"^. As in Fig. 1, three curves are shown for each model, corresponding to fractional radii 0, 1/3 and 2/3 from the cloud 
center. 



lGrvetal.Ll2002h . Figure 4 shows the effects of changing the as- 
sumed H2 formation rate constant Rg in some of the models il- 
lustrated in Fig. 1 . Although the H2 column and number densi- 
ties increase faster than linearly with changes in Rg at late times, 
the timescale may only be shortened after the molecular fraction 
approaches its maximum possible value. 

As noted in Sect. 2.4, high molecular fractions may only ac- 
crue in cases where the free-space photodissociation rate of H2 
is attenuated (by dust and self-shielding) to the point where it 
is no more than competitive with cosmic -ray induced destruc- 
tion rates. It follows that natural variations in the ambient inter- 
stellar radiation field (ISRF) may encourage or discourage H2 
accumulation, but only modestly and at late times when equi- 
librium is approached. Shown in Fig. 5 at left are results for a 
lower-density model (n(H) - 8 cm""*) at N(H) = 4 x 10^° cm"^, 
with the ISRF lowered by factors of 2 and 4. Clearly the H2 
accumulation process may be enhanced by making the material 
sufficiently dark, but only at the expense of incurring very long 
equilibrium times. Note that the darker models actually form and 
accumulate somewhat less H2 at smaller times because they are 
substantially cooler; lower photoelectic heating rates accompany 
lower photoionization rates for molecules. Clouds in very locally 
darker regions would have to be overpressured in order to com- 
pensate for this sluggishness, but in the ISM a tendency to higher 
thermal pressure at higher density seems a real effect. 



Figure 5 at right shows the effect of varying the primary cos- 
mic ray ionization rate by factors of three about the standard 
value (h - 10""' s"'. Such variations have little effect on the 
equilbrium H2 densities in standard H I clouds as long as they 
are not so large as to upset the thermal balance (which can occur 
for 10 pc). The more complicated behaviour of is 
discussed in Sect. 4. 

4. H3+ and other trace species 

4.1. Evolution of n(Hi^ ) 

In the ISM there is a general increase in N(H3^) with redden- 
ing corresponding to N(H3^) = 8 x 10'^ cm"^ over 6 magni- 
tudes, or to N(H3+)/N(H) « 2.3 x 10"*^ for the standard gas- 
redde ning ratio N(H)/Eb- v = 5.8 x 10^' cm"^ mag"' (see Fig. 
14 of iMcCall etall (|2002|) ). 

As discussed previously ( lLisztll2003l) or as shown in the fig- 
ures here, n(H3^) concentrations comparable to those required 
can be produced in equilibrium models of diffuse clouds of mod- 
erate number and column density for cosmic ray ionization rates 

at or somewhat above 10"'* s"' (see Fig. 5 at right). Such 
models of H3+ in diffuse gas have the added virtue that the 
HD/H2 ratios observed in diffuse clouds - which are very sensi- 
tive to the actual proton density - are also reproduced. 
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As shown in the figures, the H^^ fraction is sensitive to the 
same effects which increase the concentration of H2, and in some 
circumstances to the cosmic ray ionization rate as well (Fig. 5 
here). There are some more subtle effects which are less obvious 
in the figures: in models which achieve large H2-fractions, the 
electron density is lower by a factor of (typically) about two at 
late times because the gas cannot sustain high proton densities 
when most of the hydrogen is molecular. 

Although the H^^ fraction equilibrates on very much the 
same final timescale as that of H2, the problem of timescales 
is somewhat more critical for because its time derivative is 
much steeper; the final concentration occurs more nearly at only 
the final moments of the calculation. Conversely, the H^^ frac- 
tion would deteriorate more rapidly should conditions suddenly 
become less hospitable. Of course this effect is mitigated slightly 
when the fraction is calculated with respect to H2 (instead 
of n(H)) but in most cases the lines of sight with known N(H3+) 
are not accessible to direct measurement of N(H2) or even N(H 
I). Note that is also more centrally concentrated than H2, but 
the effect is really strong only in models within which there are 
larger radial gradients in n(H2). 

4.2. Effect of changing 

The presence of high cosmic ray ionization rates in difiiise gas 
has been inferred from very elementary analys is of the 
abundance which treats it like an atomic ion dMcCall et al.L 
I2OO2 I). Assuming - as is true of dark, dense molecular gas - that 
every primary ionization of H2 leads to an H3^ ion and that H3^ 
ions are destroyed by electron recombination with a rate con- 
stant a^, there follows a simple, formal expression for the 
density: 

«(H3+) * 2^„n(H2)/(a, n(e)) (3) 

The discrepancy between observation and the values of 
n(H3^) which result from putting typical values for (h, n(e), etc. 
in this expression is sufficient to suggest that » 10"'^ s"'. 
Nonetheless, Eq. 3 is only a formal, very approximate solution 
for n(H3^) because of the implicit dependences of n(e) and 
upon the physical effects of ionization and heating (respectively) 
denoted by (h, and because of other, less obvious assumptions 
made about the molecular chemistry. Furthermore, the obvious 
changes in the H3^/H2 ratio with time show that such consider- 
ations might at best apply only at very late times (e.g Fig. 5 at 
right). 

Figure 5 at right exhibits the further surprising result that the 
H3^ abundance is largely independent of the cosmic ray ioniza- 
tion rate at early times when the H2 abundance is small. This 
occurs because formation of H3^ is a multistep process with a 
resultant inverse dependence on the electron density which may 
be much faster than linear (as suggested by Eq. 3), combined 
with the tendency (suppressed but not eliminated by neutraliza- 
tion of atomic ions on small grains) for the ionization fraction to 
increase with (h- 

In a somewhat simplified form neglecting chemical reactions 
of H2^ and H3'*" with neutral atoms etc. (appropriate in diffuse 
gas), and the complexities of grain neutralization, the main steps 
involved in forming and destroying H3^ and its immediate an- 
tecedent H2^ are: 

cosmic ray + H2 — > H2^ + e, (4a) 
H2+ + e ^ 2H, iAb) 



H2+ -H H2 ^ H3+ -I- H, (4c) 

H3+ + e ^ 3H, etc (Ad) 

So, whenever the recombination described in reaction 4b 
dominates over the protonation in reaction 4c, the abundance of 
H2^ will be inversely proportional to n(e) and n(H3^) oc n(H2) 
n(H2^)/n(e) oc f/^/n(e)^. Calculating the ionization balance in an 
atomic gas according to the prescription in lSpitzeJ (Il978h indeed 
yields n(e) oc for (u> 10"'^ s"' under conditions typical of 
H I clouds. 

For the model of a typical H I cloud shown at right in 
Fig. 5, such considerations are apparently sufficient to remove 
much of the dependence on from the H3^ density at early 
times. To ascertain the conditions under which this effect is 
important, note that the recombination coefficient of H2^ is 
= 4 X 10"^/ V^K cm^ s~' and reaction 4c proceeds with a 
rate constant k - 2 08 X 10' ^ cm^ s ' in the UMIS T database, 
see iLeTeuff^etal] (l2000h or \http://w ww.rate99.co.uk Thus for 
Tk = 100 K and taking n(e) as a small multiple (2x) of the gas- 
phase abundance of carbon, n(C^) a; 1.4 x 10""* n(H), it follows 
that recombination dominates and the functional dependence of 
n(H3^) upon is weakened for n(H2)/n(H) ^ 0.05, with higher 
H2 concentrations required for higher n(e) and/or ^h- In Fig- 5 at 
right, this is case for the first 1 - 2 x lO*" yr. 

4.3. Influence of H3+ on other trace molecule abundances 

In general the surprisingly high abundance of H3^ in diffuse gas 
should not alter significantly the patterns of abundance of other 
trace species under the quiescent conditions discussed here, and 
it will certainly not by itself explain such outstanding mysteries 
as the high abundances of CH^, NH, HCO^, etc. For the former, 
the reaction C + H3^— > CH^ + H2 is exothermic but not partic- 
ularly fast and the fraction of gas-phase carbon which is neutral 
is very small in diffuse gas, typically below 1%: this forms only 
miniscule amounts of CH+ (more CH^ might be form ed by the 
reaction of C with H 2+). For NH jCrawford & Wilhamsi 1 1 9971: 
lMever&Rotail991h . the reaction of N and H3^ is endothermic. 
For HCO^, which in dark clouds is formed via the reaction H3^ 
-H CO ^ HCO+ + H2, the amounts of both H3+ and CO are too 
small in diffuse clouds to sustain HCO^ against its recombina- 
tion to CO, HCO+ + e - -^ CO -H H2. C O can form from HCO+ 
in quiescent diffuse gas (iLiszt & Lucasl I2OOQ) . but not the other 
way around. 

One species for which might make a difference is OH, 
because both the reactions O -1- H3^ — » OH^ + H2 (rate constant 
Ati = 8 X 10"'° cm 3 s"') and 0+ -H H2 ^ OH+ + H (rate con- 
stant k2 - 1.7 X 10"' cm ^ s"') are exothermic. Figure 6 shows 
the time evolution of the volume formation rate of OH^ through 
both routes (shown are fcin(0)n(H3^) and A;2n(0^)n(H2)) in three 
models threaded by a high cosmic ray flux, = 3 x 10"'^ s"', 
at three densities n(H) - yi cm"^, 64 cm~^ and 128 cm"^. 
The cloud diameter is held fixed so that N(H) = 4 x 10^*' cm"^, 
8 X 10^° cm"' and 16 x 10'" cm"', respectively, purposely driv- 
ing the highest density gas to very high H2 fractions and H3^ 
abundances at late times. 

As shown in Fig. 6, the reaction of O and H3^ is never a 
dominant source of OH^. The figure does show, however, the 
effect of io nizatio n balance on the OH formation rate (see also 
Sect. 2.3 of lLiszj (12003)): at later times the first step toward OH 
formation occurs more vigorously in the more tenuous models 
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time [yr] time [yr] 



Fig. 4. Similar to Fig. 1, for models having n(H) - 32 cm ^, N(H) = 4 X 10^° cm ^ and primary cosmic ray ionization rate 
= 10"'^ s"', with the rate constant for H2 formation scaled by factors of 2 about the standard value. 




10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 
time [yr] time [yr] 

Fig. 5. H2 and Ha^ densities calculated for varying strengths of the interstellar radiation field (left) and cosmic ray primary ionization 
rate (right): the scale for n(H2) is given at far left, that for n(H3^) at far right. In the panel at left the scaling parameter Go is 
taken as 1, 1/2 and 1/4; at right varies by factors of 3, 1 and 1/3 about the standard model value = 10"'^ s"'. As in Fig. 1, 
three curves are shown for each model, at fractional radii 0, 1/3, and 2/3 from the center. The models in both panels have N(H) 
= 4 X 10^'^ cm"^ but the model at left has lower number density n(H) = 8 cm"^. Note that n(H3^) is relatively insensitive to changes 
in at early times when n(H2) is small, as discussed in Sect. 4 of the text. 
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10* 10^ 10^ 10^ 
time [yr] 



Fig. 6. OH^ volume formation rates due to reaction of + H2 
(upper family of curves) and O + Hj^ (lower family) in the pres- 
ence of a high primary cosmic ray ionization rate = 3 x 10""' 
s Results are shown at the usual three fractional radii within 
models of diameter 4.05 pc having n(H) - 32 cm"-', 64 cm"-' 
and 128 cm^^ , so that N(H) = 4 x 10^" cm^^, 8 x 10-" cm^^ 
and 16 x 10^" cm"-, respectively). 



because their higher ionization levels (both overall and in oxy- 
gen) more than compensate for lower H2 fractions. In the mod- 
els with larger H2 fractions, the proton density declines at late 
times because a more nearly molecular gas cannot sustain a high 
proton density and the 0"^/0 ratio also declines, due partly to its 
rapid charge exchange reaction with hydrogen and because of its 
rapid reaction with H2. The cosmic-ray ionization alone, even at 
the high rates assumed here, cannot support high ionization frac- 
tions in atomic hydrogen or oxygen in the presence of PAH and 
high H2 concentrations. 

The next (necessary) step toward OH formation, OH+ -H 
H2— » OH2^ + H will dominate over recombination of the OH^ 
ion for H2-fractions as small as 2% at diffuse cloud temperatures. 
Thus, the thinner models will actually form OH more rapidly 
(not just OH^), as the H2 accumulation process equilibrates 0. 
Chemical effects like those shown in Fig. 6 undoubtedly con- 
tribute to the appearance of radio frequency OH emission in the 
envelopes of d ark clouds well outsid e the most st rongly molecu- 
lar reg ions (see lWannier et alj(ll993h and Fig. 6 ofLiszt & LucasI 
(Il996l) ). In the current context we expect/predict that regions of 
higher than normal OH abundance should be found optically in 
diffuse clouds having modest extinction and H2fractions. 



^ At early times, the formation of OH in less dense clouds -will be 
seriously impeded by their larger e/H2 ratios and the denser models may 
make more OH. As with H3+ itself, high OH abundances only appear at 
the conclusion of the H2 accumulation process 



5. Summary and discussion 

5. 1 . Summary 

We showed the approach to equilibri um under benign conditions 
for some recently-published models ( Liszt & LucasLl2000tlLisztL 
|2()02, 2003), calculating the time-dependent H2 and Hs'*^ forma- 
tion in pc-size, otherwise largely static gas parcels meant to rep- 
resent typical quiescent diffuse "clouds". We pointed out that 
self-shielding does not hasten the formation or accumulation of 
H2, although it may foster the production of more H2 over longer 
times. Because the H2 formation process in diffuse gas is incom- 
plete, timescales can in general not be deduced from simple scal- 
ing arguments based on the local H2 formation rate. The equilib- 
rium H2 formation timescale in an H I cloud having a column 
density N(H) a; 4 x 10^° cm"^ is 1 - 3 x 10^ yr, nearly inde- 
pendent of the assumed density or H2 formation rate on grains, 
etc. Under the benign conditions described here, the bulk of the 
H2 formation in the diffuse ISM actually occurs under diffuse 
conditions in physically larger regions of moderate density; that 
is, there is little apparent advantage in hypothesizing a process 
which forms H2 in denser clumps, if the density remains within 
the range inferred for quiescent diffuse gas. 

Other species reach equilibrium essentially as soon as H2 
does, even when they involve higher-order and/or very "slow" 
processes (like cosmic ray ionization of H2) but their equilib- 
rium abundances appear more nearly at just the last minute; it is 
more critical for explaining Ha^ that conditions equilibrate than, 
say, for H2 itself. This is also relevant because the chemistry of 
a weakly-molecular gas - which could be either a model of low 
density or one of higher density at early time - behaves in var- 
ious particular ways. For instance, the predicted abundance of 
H3^ is less strongly dependent on the cosmic -ray ionization rate 
when n(H2)/n(H) ^ 0.05. 

Attempts to form increased amounts of H3+ may be thwarted 
at very high assumed values of both because the fractional 
ionization increases and because the H2-fraction declines as 
cosmic-ray ionization overtakes photodissociation as the main 
H2 -destruction mechanism. Assuming (h» 10"'^s"' will also 
change the thermal balance in the diffuse ISM as a whole as 
cosmic-ray ionization becomes the dominant heating process 
in both warm and cold neutral diffuse gas. Although observa- 
tions of both HD and H3^ seem to indicate an enhanced low- 
level source of hydrogen ionization, it is important to ascer- 
tain whether it is indeed a cosmic-ray related process which 
is responsible (or some other form of very hard radiation) and 
whether any related effects occur in other phases of the interstel- 
lar medium, either in the warm or dark gas. 

The presence of relatively high H^^ abundances per se prob- 
ably does not have strong effects on the chemistry of other trace 
species whose abundances are puzzling (CH+, NH, HCO^ etc.) 
or on those of species like OH which are generally considered to 
be better understood. However, in considering competition be- 
tween OH formation paths through H2 and Ha^, we showed that 
late-time OH formation proceeds most vigourously in more dif- 
fuse regions having modest density, extinction and H2 fraction 
and somewhat higher fractional ionization, suggesting that high 
OH abundances might be found optically in diffuse clouds hav- 
ing modest extinction. 

5.2. Lifetimes of diffuse clouds and their H2 

In this work we modelled the formation and accumulation of 
H2 under conditions like those which are inferred from observa- 
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tions of the H2 and related neutral gas species themselves; tem- 
peratures from H2 dSavage et all 119771: |Shull et al.', '2004") and 
meari thermal pressures (densities) from C I ({Jenkins & Tripp, 
l200lh . Although modest amounts of H2 form quickly enough 
that molecule-free sightlines are relatively rare, the observed 
high mean molecular fractions of 25%-40% jSavage et al.Lll977t 
iLiszt & Lucasll200 2') are achieved only over timescales of 10^ yr 
or longer, when starting from scratch. 

H2 aside, timescales for the r mal, i onization, chemical 
etc. equihbrium are s hort dSpitzed 119781: IWolfire et all 120031: 
iRawhngs et al.L |2002|) and diffuse gas is not normally thought 
of as having much longevity in any particular form. The H I 
cloud spectrum evolves under the continual influences of cloud- 
cloud collisions, evaporation , passage of supernova bla st waves 
and stellar wind shocks etc. dChieze & Lazaret 1198 0') and in a 
3-phase model of the diffuse ISM, gas in any particular volume 
el ement turns over on times cales shorter than 10'' yr according 
to lMcKee & Ostriked(ll977l) . 

How then are high H2 fractions and short diffuse gas life- 
times to be reconciled? This question has never really been 
considered when the overall structure of the diffuse ISM is 
modelled, but there are several alternatives. 1) The long equi- 
librium timescales under diffuse neutral conditions provide, in 
reverse, some resilience which might help pre-existing H2 to 
weather the changes which occur 2) Dark clouds have haloes 
(Wannieret al., 1993; Bensch, 2006) and it has been suggested 
that the molecules supposedly seen in diffuse clouds are actu- 
ally located around and e xchanging material with dark clouds 
jFederman & Allenl Il99lh : this may be true in some cases, but 
complex trace molecules are also seen along sightlin es which are 
quite well separated from the nearest dark material (i Liszt et all 
|2002|) . 3) Perhaps most promisingly, turnover within diffuse neu- 
tral gas due to turbulence, rapidly cycling material through a 
very dense phase, ma y provide for much faste r formation of H2 
and t race species jGlover & Mac Lowl 120051: 'Falgaron e et al.L 
l2005h . establishing high molecular fractions from scratch on 
short timescales. The small fraction of neutral material which 
is observed to exist at very high pressure within diffuse clouds 
jjenkins & Tr ipp, 20 1]) and other asp ects of small-scale struc- 
ture in diffuse gas (IDeshpandet l2000l) may perhaps also be un- 
derstood in such terms. 
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